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Abstract 

The long standing problem of non perturbative renormalization of a gauge field 
theoretical Hamiltonian is addressed and explicitly carried out within an (effec- 
tive) light-cone Hamiltonian approach to QCD. The procedure is in line with 
the conventional ideas: The Hamiltonian is first regulated by suitable cut-off 
functions, and subsequently renormalized by suitable counter terms to make it 
cut-off independent. The formalism is applied to physical mesons with a dif- 
ferent flavor of quark and anti-quark. The excitation spectrum of the p-meson 
with its excellent agreement between theory and experiment is discussed as a 
pedagogical example. 

1. Introduction 

When starting in 1984 with Discretized Light-Cone Quantization (DLCQ) [ 1 1 and with a revival of 
Dirac's Hamiltonian front form dynamics [2], all challenges of a gauge field Hamiltonian theory were 
essentially open questions, particularly the non perturbative bound state problem, the many -body aspects, 
regularization, renormalization, confinement, chirality, vacuum structure and condensates, just to name 
a few. The step from the gauge field QCD Lagrangian down to a non relativistic Schrodinger equation 
was completely mysterious. Now we know better [3 ]. We have understood, for example, that the chiral 
phase transition, in which the quarks are supposed to get their mass, is not the major challenge. The 
challenge is to understand what happens after the phase transition, at zero temperature. The challenge 
is to understand the spectrum of physical hadrons and to get the corresponding eigenfunctions, the light 
cone wave functions. The light-cone wave functions \£ for a hadron with mass M encode all possible 
quark and gluon momentum, helicity and flavor correlations and, in principle, are obtained by diagonal- 
izing the QCD light-cone Hamiltonian H LC \^i) = M?|^), H LC = P + P~ - P±, in a complete basis 
of Fock states with increasing complexity. For example, the positive pion has the Fock expansion: 

|^+) = $>|7r+)|n) = k ±l )\ud) + * { u % /7r (x t , k ± i)\udg) + ... , 

n 

representing the expansion of the exact QCD eigenstate at scale A in terms of non-interacting quarks and 
gluons. The particles in a Fock state have longitudinal light-cone momentum fractions xi = kf / P + 
and relative transverse momenta k^. The form of ^ n / H (xi, k±i) is invariant under longitudinal and 
transverse boosts; i.e., the light-cone wave functions expressed in the relative coordinates x-i and k±i are 
independent of the total momentum (P + , P±) of the hadron. The first term in the expansion is referred 
to as the valence Fock state, as it relates to the hadronic description in the constituent quark model. The 
higher terms are related to the sea components of the hadronic structure. It has been shown that the rest 
of the light-cone wave function is determined once the valence Fock state is known with explicit 

expressions given in [5]. 

The key issue is to overcome the problem of any gauge theory, that the unregulated theory exposes 
logarithmic singularities. The problem of regularization and renormalization has been solved in the 
perturbative context of scattering theory, but not in the non perturbative context of a Hamiltonian. It is 
addressed to in the first two sections and applied in the remainder of this paper. 




Fig. 1 : Regularization of the interaction by vertex regularization. In a matrix element, as illustrated on the left for a vertex, a 
quark changes its four-momentum from fei to fea, thus Q 2 — — (Jfej — fe) 2 • The vertex interaction is regulated by multiplying 
with a form factor F(Q 2 ), as indicated by the circle. — Instantaneous interactions are treated correspondingly, as shown on 
the right for a seagull. 

2. Regularization 

Canonical field theory with the conventional QCD Lagrangian allows to derive the components of the 
total canonical four-momentum P 11 . Its front form version [3| rests on two assumptions, the light cone 
gauge A + = [ 6 1 and the suppression of all zero modes (3J0- The front form vacuum is then trivial. 
I find it helpful to discuss the problem in terms of DLCQ In the back of my mind I visualize 

an explicit finite dimensional matrix representation of the Light-Cone Hamiltonian as it occurs for finite 
harmonic resolution. Such one is schematically displayed in Fig. 2 of 0. All of its matrix elements are 
finite for any finite x and k±. 

The problem arises for ever increasing harmonic resolution, on the way to the continuum limit: 
The numerical eigenvalues are numerically unstable and diverge logarithmically fT6l . contrary to the 
calculations in 1+1 dimension One must regulate the theory. 

As usual, regularization is not unique, and can be done in many ways. But not all regularization 
techniques of the past are applicable in an Hamiltonian approach. Dimensional regularization, for ex- 
ample, is not applicable in a matrix approach which is stuck with the precisely 3+ 1 dimensions of the 
physical world. Perturbative regularization [ 8 1 is not applicable in the non-perturbative context. The 
Fock space regularization of Lepage and Brodsky [ 6 1 , see also [ 3 1 has blocked renormalization for many 
years. It was impossible to find suitable counter terms. In recent work, the invariant mass squared regu- 
larization has been abandoned in favor of Pauli-Villars regularization [9 10 1. But thus far it is unclear 
how the approach 1 10 1 is applicable to the spectra of physical mesons. After applauding the light-cone 
approach 1111 . Wilson and collaborators lH2l have attempted to base their considerations almost entirely 
on a renormalization group analysis, but no concrete technology has emerged thus far. 

In recent years I have favored an other regularization, which allows for renormalization and which 
is technically sufficiently simple to be carried out explicitly. All divergences can be traced back to 
Dirac's relativistic vertex interaction (k\, hi\V\k2, h 2 \ k%, h 3 ), in which some particle '1' is scattered 
into two particles '2' and '3' with their respective four-momenta k and helicities h, see Fig. [2 The 
matrix element for bremsstrahlung, for example, is proportional to k±, (&i,T |^|&2> T> &i> T) °t |^l[> 
see Table 9 in [3], when the quark maintains its helicity while irradiating a gluon with four-momentum 
fcg = (xP + , k±,kz). Singularities arise typically when squares of such matrix elements are integrated 
over all k± as in the integrations of perturbation theory. 

The singularities are avoided a priori by vertex regularization, by multiplying each (typically off- 
diagonal) matrix element with a regulating form factor F: 

(k 1 ,h 1 \V\k 2 ,h 2 ;k 3 ,h 3 ) (k 1 ,hx\V\k 2 ,h 2 ;h,h 3 ) F(Q) . (1) 

It took several years to realize that it is the Feynman four-momentum transfer across a vertex, Q 2 = 
— (fci — k 2 ) , which governs any effective interaction. The minimal requirement for such a form factor 
is limQ_>o F{Q] A) = 1 and limq^oo F(Q; A) = 0. The job would be done by a step function, F(Q) = 
Q(Q 2 — A 2 ). The limit A — > suppresses the interaction all together, the limit A — > oo restores 



the interaction and its problems. Any finite value of A 2 restricts Q 2 to be finite and eliminates the 
singularities. But the sharp cut-off generates problems in an other corner of the theory and F(Q) must 
be an analytic function of Q, as to be seen below. 

Vertex regularization takes thus care of the ultraviolet divergences. The (light-cone) infrared sin- 
gularities are taken care of as usual by a kinematical gluon mass. 



3. Renormalization 

The non perturbative renormalization of the Hamiltonian was stuck for many years by the fact that the 
coupling constant g and the regulator function F(Q) multiply each other in Eq.Q. It was always clear 
that one may add non local counter terms fT2l . but is was not clear how they could be constructed. 
Progress has come from recent work on a particular model [19], which did allow to formulate a paradig- 
matic example in modern renormalization theory. 

Here is the general but abstract procedure. 

Suppose to have solved i^Lcl^i) = Mf\^>i) for a fixed value of the 7 'bare' parameters in the 
Lagrangian, for the coupling constant g = go and the 6 flavor quark masses rrif = (m/)o, and for a 
fixed value of exterior cut-off scale A = Ao- Suppose further that these 7+1 parameters are chosen such, 
that the calculated Mf agree with the corresponding experimental values. Next, suppose to change the 
cut-off by a small amount 5 A. Every calculated eigenvalue will then change by <5M 2 . Renormalization 
theory is then the attempt to reformulate the Hamiltonian, such, that all changes 5Mf vanish identically. 
The fundamental renormalization group equation is therefore for all eigenstates i: 



dM„ 



1 o 



dMi 



= . (2) 

g=g ,m f =m fo ,A=A 



Equivalently one requires that the Hamiltonian is stationary, 5H^c = 0, with respect to small 5 A. 
Hence forward reference to the renormalization point (go,mf , Aq) will be suppressed. 

The Hamiltonian can be made stationary by making g and the m f functions of A, by introducing 
physical coupling constants and masses, g and mj, respectively, which themselves are functions of the 
bare g and mj, and which are functional of the regulator F = F. The variation of H^c reads then 
SHlc = 5g s + J2f d™f S 8m^ + ^i 7 * 5 ^ = 0, with the familiar variational derivatives. However, 

since g and mj are themselves functionals of F, this reduces to 5H^c = 5F SI ^jf' = 0. The fundamental 
equation of renormalization theory $2% is then replaced by 

— dF 

6F = 6A dX = °- ^ 

It can be solved by counter term technology, as follows. A counter term is added to the Hamiltonian, 
whose interaction has exactly the same structure except that the regulator F(Q) is replaced by C(Q). 
This defines 

F(Q, A) = F(Q, A) + C{Q, A) , (4) 

subject to the constraint that the counter term vanishes at the renormalization point, C(Q, A)| a=Aq = 0. 
The fundamental equation (|3} defines then a differential equation dC{Q]A)/dA = — dF(Q; A)/dA, 
which, in its integral form, includes the initial condition 

A 

C(Q, A) = - Jds dF f s ,s) = F{Q, A ) - F(Q, A) . (5) 

Ao 

The renormalized regulator function, F = F + C, is manifestly independent of A: 



F(Q, A) =F(Q, A ) 



(6) 



By construction, the value of Ao is determined by experiment. One should emphasize an important point: 
In deriving Eq.©, use was made of assuming the regulator function has well defined derivatives with 
respect to A. The theta function of the sharp cut-off, however, is a distribution with only ill defined 
derivatives. This raises an other important point: If F(Q,A) is a function of Q/A other than a theta 
function, one must specify how the function approaches the limiting values of 1 and of 0. The case of the 
'soft' regulator F(Q,A) = A 2 /(A 2 + Q 2 ) is only a very special example. In a more general approach 
the soft regulator plays the role of a generating function 



F(Q,A) 



N 



l + £(-l) n s n A n 



n=l 



gn 
9A" 



A 2 



A 2 + Q 2 



(7) 



The partials A n d n /dA n are dimensionless and independent of a change in A. The arbitrarily many 
coefficients s n are renormalization group invariants and, as such, subject to be determined by experiment. 



4. The effective (light-cone) Hamiltonian 

In a field theory, one is confronted with a many-body problem of the worst kind: Not even the particle 
number is conserved. For to formulate effective Hamiltonians more systematically, a novel many-body 
technique had to be developed, the method of iterated resolvents (5][HI> whose details are not important 
here. Important is that the effective light-cone Hamiltonian -ff e LC has the same eigenvalue as the full 
light-cone Hamiltonian H^q and that it generates the bound state wave function of valence quarks by an 
one-body integral equation in (x, k±): 



M 2 4> hlh2 (x,k±) 



m\ + k 2 m| + k 2 



x 



1 



dx'd 2 k' L *l> h > ih i 2 (x',k[_) 
y/x(l - x)x'{l - x>) 



F(Q q )F(Q 




[u(ki,hi)^u(k[,h 2 )][v(k 2 ,h' 2 )^ v (k 2 ,h 2 )] ■ (») 



One has achieved H c i,c\^qq) = M 2 \^ q q). Here, M 2 is the eigenvalue of the invariant-mass squared. 
The associated eigenf unction 7jJh 1 h 2 { x i k±) is the probability amplitude (x, k±, hi; 1 — x, — k±, h 2 \^ q q) 
for finding the quark with momentum fraction x, transversal momentum kj_ and helicity hi, and corre- 
spondingly the anti-quark. Expressions for the (effective) quark masses mi and m 2 and the (effective) 
coupling function a(Q) are given in [ 14]. Q q and Qq are the Feynman momentum transfers of quark 
and anti-quark, respectively, and u{ki, hi) and v(k 2 , h%) are their Dirac spinors in Lepage Brodsky con- 
vention [6], given explicitly in 0. Finally, the form factors F(Q) restrict the range of integration and 
regulate the interaction. Note that the equation is fully relativistic and covariant. Note also that Eq.® 
is valid only for quark and anti-quark having different flavors 004). The additional annihilation term 
for identical flavors is omitted here, and presently investigated by Krahl [15|. The 'mean momentum 
transfer', Q 2 = \ {Q 2 q + Q 2 ^, allows to replace Eq.® by 



M 2 i(j hlh2 (x,k±) 



m 2 + k 2 ^rri2 + fe£ 



1 — x 



1 rdx'd 2 k' ± VV/if { x \ k\) 

^h 1 h 2 (x,k±) -rj L / n 7^ 
i7T h^t' J ^fx(l-x)x'{l-x') 

"■q>' L q 



Q 2 



R(Q) [nik^hiyf^X)] [v(k' 2 ,h' 2 ) lfi v(k 2 ,h 2 )] 



(9) 



The form factors F{Q) have made their way into the regulator function R(Q) = F 2 {Q). Krautgartner 
and Trittmann [16] have shown how to solve numerically such an equation with a high precision. But 
since the numerical effort is so considerable, it is reasonable to work first with (over-)simplified models, 
as specified next. 



The Singlet-Triplet model. Quarks are at relative rest when kj_ = and x = mi/ (mi + W12). 
An inspection of Eq.(33) in ifOl reveals that for very small deviations from the equilibrium values, 
the spinor matrix (hi, h2\S\h' 1 ,h' 2 } = \u(ki, hi)j^u(k[, h' 2 )\ \v(k 2 , h' 2 )^^v(k 2 , h 2 )\ is proportional to 
the unit matrix, (hi,h 2 \S\h' 1 h' 2 ) ^ 4mim2 5 hl fl > 5 h2jh ^. For very large deviations, particularly for 

k^ > kl, holds Q 2 ~ k\_ 2 and (U |5| U) ^ 2fc| 2 . 
The Singlet-Triplet (ST) model combines these aspects: 



</ti,/t 2 |5|/n,fe 2 ) J ^3^ + 2, for/n = -/i 2 

Q2 i fat h x = h 



2- 



(10) 



For anti parallel helicities /ii = — hi (singlets) the model interpolates between two extremes: For small 
momentum transfer Q, the '2' in Eq. (fT0l is unimportant and the Coulomb aspects of the first term prevail. 
For large Q, the Coulomb aspects are unimportant and the hyperfine interaction is dominant. For parallel 
helicities hi = h 2 (triplets) the model reduces to the Coulomb kernel. The model over emphasizes many 
aspects but its simplicity has proven useful for fast and analytical calculations. Most importantly, the 
model allows to drop the helicity summations which technically simplifies the problem enormously. The 
model can not be justified in the sense of an approximation, but it emphasizes the point that the '2', or 
any other constant in the kernel of an integral equation, leads to numerically undefined equations and 
thus singularities. Replacing the function a(Q) by the strong coupling constant a s = g 2 /Air completes 
the model assumptions. Hence forward, the overline bars for the effective quantities will be suppressed. 

5. The potential energy 

It is possible to subtract a c-number from H e ^c and to define an effective Hamiltonian H c g implicitly by 

H cLC = (mi + m 2 ) 2 + 2 (mi + m 2 ) H eS , H eS \<p) = E\tp) . (11) 

Its eigenvalues have the dimension of an energy. Note that mass and energy in the front form, on the 
light cone, are related by 

M 2 = (mi + m 2 ) 2 + 2 (mi + m 2 ) E , (12) 

and not by M 2 = (mi + m,2) 2 + 2 (mi + m 2 ) E + E 2 , as usual. Only if the energy is negligible as 
compared to the quark masses, i.e. only if (E/ (mi + m 2 )) 2 1, the two relations coincide. 

A rather drastic technical simplification is achieved by a transformation of the integration vari- 
able. One can substitute the integration variable x by the integration variable k z , which, for all practical 
purposes, can be interpreted [3| as the z-component of a 3-momentum vector p = (k z , k±). For equal 
masses mi = m 2 = m, the transformation is, together with its inverse, 



1 + 



(x-±) 2 

k 2 (x) = (m 2 + k 2 ) \ 2J . (13) 
x{l — X) 



'm 2 + k 2 + k 2 

Inserting these substitutions into Eq.© and defining the reduced wave function ip by 




jA(k z ,k±) _ 

^h!h 2 (x,k±) = Vh x h 2 (k z ,ki_) , A(p) = \l-\ 2' ( 14 ) 

^£(1 — x) V m 

leads to an integral equation in the components of p, in which all reference to light-cone variables has 
disappeared. Using in addition the ST-model of Eq.lfTOT). Eq.© translates for singlets identically into 



d 3 p' ( Am 2 \ R(Q) 



M 2 cp(p) = 4 m 2 +p 2 <p(p) - -\ \ ^=4== — ^ + 2 ip(p') , (15) 



2vr 2 J y/A(p)A(j/) \ Q 



m 



f a c f d 3 p' /4m 2 \ R(Q) 



with a c = |a s . The equation for the triplets is obtained by dropping the '2'. In the ST-model, the 
helicity arguments in the wave functions can be suppressed. Applying the relation between mass and 
energy, as given in Eq. (fT2l . the equation is converted to 

f ,^ a c f d 3 p' ( 4m 2 ( o \ R(Q) 

since the reduced mass for mi = mi = m is m r = m/2. The first term in this equation, p 2 /2m r , 
coincides with the kinetic energy in a conventional non-relativistic Hamiltonian. This is remarkable in 
view of the fact that no approximation to this extent has been made. The fully relativistic and covariant 
light-cone approach has no relativistic corrections in the kinetic energy! 

Since the first term in Eq. dToT i is a kinetic energy, the second must be a potential energy — in a 
momentum representation. In principle, it could be Fourier transformed with e~ ipr to a configuration 
space with the variable r. But due to the factor A{p)A{p') in the kernel, the resulting potential energy 
would be non-local. The non-locality of the potential is certainly mathematically exact. But I do not 
expect this to generate aspects of leading importance, and avoid it by the simplification A{p) = 1, 
both in Eqs. (fT4b and (fTBT) . With A(p) = 1, the mean four momentum transfer Q 2 reduces to the three 
momentum transfer q 2 = (jp — p') 2 . In consequence, the kernel of Eq. lfToT) depends only on q = p — p', 

+ nr>/A.^«, (17) 

Its Fourier transform is a local function, which plays the role of a conventional potential energy V(r) in 
the Fourier transform of Eq. JToT l. i.e. in 



E ip(r) 



p 2 



i>(f) ■ (18) 



Here is the simple Schrodinger equation ! Despite its conventional structure it is a front form equation, 
designed to calculate the light-cone wave function ip(f) — > ip(p) — > ip q q{x, k±). I conclude this section 
with a subtle point, which needs clarification in the future. The simplification A(p) = I is different from 
a non-relativistic approximation. The approach is certainly valid also for relativistic momenta p 2 3> m 2 , 
particularly Eqs. dToT i and fffl . The reason is that A(jp) occurs only under the integral. There, the large 
momenta are suppressed by the regulator, anyway. 



6. The renormalized Coulomb potential 

Hence forward, I restrict consideration to the triplet case, i.e. to Coulomb kernels like U(q) ~ R(q)/q 2 . 
The point is that he renormalized Coulomb potential is always finite at the origin, as opposed to the 
conventional Coulomb potential with its ^-singularity. It is instructive to verify this explicitly for the 
sharp cut-off as a regulator, that is for U(q) = — 2 ^ c q i @{q 2 — ^ 2 )- The Fourier transform according 
to Eq. dTTl gives V(r) = — ^ -Si(Ar). Using the well known asymptotic expansions of the Integral 

Sine Si(Ar) gives lim r ^oo V(r) = -2c and lim r _* V(r) = a c A(-| + ^^)- The regulated Coulomb 
potential is finite at the origin. The cut-off dependence near the origin is one of the most important 
insights of the present work and has a deep physical reason to be discussed below. In analogy to Eq.Q 
the regulator is chosen as 



R(q) 



N 



l + J2(-l) n S n X r 



n=l 



Qn 



A 



\ 2 +q 2 



which gives straightforwardly the generalized Coulomb potential 



V(r) 



N 



n=l 



1 + e 



-Ar 



n=0 



(19) 
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Fig. 2: Schematic behavior of the renormalized Coulomb po- 
tential, see also the discussion in the text. 
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Fig. 3: The dimensionless Coulomb potential Wn(jj\ 0, 1, 0) 
is plotted versus the radius parameter y — Xr for different N, 
i.e. from bottom to top for TV = 4, 5, 6, 7, 8. 



with so = 1. This result illustrates an other important point: The Laguerre polynomials are a complete set 
of functions. The term added to the -1 in Eq . d20t is thus potentially able to reproduce an arbitrary function 
of r. The description in terms of a generating function is therefore complete. It is often useful to study 
V(r) in the dimensionless form by introducing Wn{u] {s n }) = V(r)/a c X = ( — l+e~ y Y,n=o s nV n )/y- 
It is a function of r only through the dimensionless y = Xr and depends parametrically on the N 
coefficients {s n }. 

The physical picture which develops is illustrated in Fig. |2] In the far zone, for sufficiently large 
r, the potential energy coincides with the conventional Coulomb potential — — . Since the potential is 
attractive, it can host bound states which are probably those realized in weak binding. In the near zone, 
for sufficiently small r, the potential behaves like a power series CQ+c\r+C2r 2 which potentially can host 
the bound states of strong coupling. Since Eq.(EUt is an analytic function of r, the actual potential must 
interpolate between these two extremes in an intermediate zone, for example by developing a barrier 
of finite height. The onset of the near and intermediate regimes must occur for relative distances of the 
quarks, which are comparable to the Compton wave length associated with their reduced mass. If the 
distance is smaller, one expects deviations from the classical regime by elementary considerations on 
quantum mechanics, indeed. 

The large number of parameter in Eq.(l2Ul can be controlled by expressing all coefficients {s n } 
in terms of three parameters a, b, and c: 

1 ^ a b c 

Sn ~ri + (n-1)! + (n - 2)! + (n-3)! ' ( } 

The first 3 coefficients are then explicitly so = 1, si = 1 + a, and S2 = \ + a + b. The dimensionless 
Coulomb potential depends then only on three parameters: Wn(v', a, 6, c). In the near zone, it is at most 
a quadratic function of y, Wiq{y; a,b,c) = a + by + cy 2 , independent of N. The remainder starts at 
most with power y N+1 . A value of a = c = and 6=1 should therefore yield a linear set of functions 
Wn(v', b, c) = y in the near zone. As shown in Fig.|3]this happens to be true for surprisingly large 
values of y. The value of N essentially controls the height of the barrier. Similarly, Wjv(y; 0, 0, 1) = y 2 
generates a set of functions which are strictly quadratic in the near zone. Again, controls the height of 
the barrier, as seen in Fig. |5] 

7. Determining the parameters by experiment 

The QCD-inspired model developed thus far has a considerable number of renormalization group in- 
variant parameters, which must be determined once and for all by experiment. In doing this [2Q|, we 




1 2 3 4 5 6 



n 

Fig. 4: The invariant mass-squares of all available tv + - 
and p + -states are plotted versus a counting index n. The 
straight lines correspond to M% = Mq + n%, with the value 
X = 1-39 GeV 2 , taken from Anisovich et al. 1228 . The filled 
circles correspond to states which have been seen empirically 
1211 . the empty ones correspond to the predictions 1221 . — 
Plot courtesy of Shan-Gui Zhou. 
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Fig. 5: The continuous lines display the generalized Cou- 
lomb potential V(r) = Q c AWjv(Ar; a, 0, c) in physical units 
as function of r, for the values N = 4, 5, 6 from bottom to 
top. The dashed line displays the harmonic approximation. 
The horizontal lines on the left indicate the oscillator states. — 
The circles indicate the experimental eigenvalues E n for the 
p + . They agree with the calculated eigenvalues for N = 6, 
shown by the horizontal lines. See the discussion in the text. 



have been inspired by the work of Anisovich et al. [22]. Enumerating the excited states of a hadron 
by a counting index n = 0,1,2,..., these authors have found the linear relation M% = Mq + n\ for 
practically all hadrons. As an example, I present in Fig.|4]the spectrum of the it- and the p-meson. The 
linear relation between mass-squared and energy on the light cone, Eq. (fT2l . allows then to conclude that 
the potential energy in the near zone must be a pure oscillator, 

V(r) = -c t + ^ftr 2 , (22) 

at least to first approximation, thus b = 0. If one addresses to reproduce the spectra of all flavor off- 
diagonal triplet mesons (pseudo-vector mesons), except the topped ones, one has to determine 6 pa- 
rameters: The 2 constants from the oscillator model, q and f t , and the 4 effective flavor quark masses 
m u = md, m s , m c , and m&. To determine them, one needs 6 experimental triplet masses. I take from 
I2T1 : The ground and first excited states of the ud and us mesons, and the ground states of uc and ub. 
The so obtained parameter values are: 

m u = 0.218 GeV, m s = 0.438 GeV, m c = 1.749 GeV, m b = 5.068 GeV, 
c t = 0.880 GeV, f t = 0.0869 GeV 3 . ( } 

The numbers differ slightly from those in [20], due to choosing a different set of empirical data, but they 
yield about the same overall agreement with all available experimental states of pseudo-vector mesons. 

Reverting the argument, one concludes as in [20] that the oscillator model in Eq. d22t explains 
quite naturally the systematics found by Anisovich et al. [22|. 

8. Relating the oscillator model to QCD 

The oscillator model in Eq. d22t is only the harmonic approximation to the QCD-inspired, generalized 
Coulomb potential in Eq. d20l . Their parameters are related obviously by q = — a c Xa, 6 = 0, and 



ft = 2a c \ 3 c. One needs more experimental information to pin down the value of a, c and N. Choosing A 
as the QCD scale, i.e. A = 200 MeV, one can use the expressions for a(Q) in [ 14] to calculate a s = a(0) 
from the measured value of the coupling constant at the Z-mass Mz = 91.2 GeV, a(Mz) = 0.118, thus 
a s = a(0) = 0.1695. Having fixed and A allows to calculate a and c from q and f t , i.e. 

a = —19.5 and c = 24.0, one can draw the generalized Coulomb potential V(r) = a c XWjsf(Xr; a, 0, c) 
for different N as done in Fig. |5] The 'experimental' eigenvalues Eq,E%,E\, and E$ for the p-meson 
are obtained from E n = (M% — 4m^)/(4m u ) and also inserted together with the empirical limits of 
error. The experimental error 5E p ^ ~ ±0.5 GeV (thus 5M Pt s ~ ±0.1 GeV) is hypothetical, since M p ^ 
is not confirmed. Taking it for granted, the lowest possible value for N is thus N = 6. This completes 
the determination of all parameters. They are universal within the model. I thank Omer[23 ] for the exact 
eigenvalues prior to publication. 

9. Summary and Conclusions 

This work is an important mile stone on the long way from the canonical Lagrangian for quantum chromo 
dynamics down to the composition of physical hadrons in terms of their constituting quarks and gluons, 
by the eigenf unctions of a Hamiltonian. As part of a on-going effort [24|, a denumerable number of 
simplifying assumptions had to be phrased for getting a manageable formalism. Among them is the 
formulation of an effective interaction by the method of iterated resolvents. As long as the assumption 
are not proven at least a posteriori, one must speak of an approach inspired by QCD. The probably 
strongest assumption in the present work, the simplifying Singlet-Triplet model, is also the one which 
can be relaxed the easiest in upcoming work. But its simplicity has been an advantage to unreveal the 
physical content of gauge theory by analytical relations. 

The biggest progress of the present work is related to a consistent regularization and renormaliza- 
tion of a gauge theory. The ultraviolet divergences in gauge theory are caused less by the possibly large 
momenta of the constituent particles, but by the large momentum transfers in the interaction. In a Hamil- 
tonian approach, such as the present, one has not much choice else than to chop them off by a regulating 
form factor in the elementary vertex interaction. The form factor makes its way into a regulator function 
which suppresses the large momentum transfers in the Fourier transform of the Coulomb interaction. 

The arbitrariness in chopping off the large momentum transfers is reflected in the arbitrariness of 
the potential at small relative distances. It is this arbitrariness which allows for a potential pocket which 
binds the quarks in a hadron. The problem how to fix this arbitrariness by experiment, in practice, is less 
difficult than anti-cipated. It suffices to determine only three parameters: a, c, and N. 

The potential energy of the present work vanishes at an infinite separation of the quarks. This 
seems be be in conflict with the potential energies of phenomenological models which rise forever. It 
also seems to be in conflict with lattice gauge calculations. Is a finite ionization limit in conflict also 
with 'confinement', i.e. with the empirical fact that free quarks have not been observed? — The present 
model prohibits free quarks as a stable solution, since the sum of the constituent quark masses is always 
larger than the mass of the corresponding hadron and a pion. Free constituent quarks would hadronize 
very quickly into bound states. This is different from atomic physics with its free constituents, where the 
binding energy is always much smaller than the mass of positronium proper. 

The most disturbing aspect of the present work is its obvious conflict with lattice gauge calcula- 
tions. I have not checked to which extent a possible linear term in the potential spoils the present excellent 
agreement between theory and experiment. The calculation of the potential energy on the lattice rests on 
the assumptions of static quarks, of quarks with an infinitely large mass. Even with present day comput- 
ers lattice gauge calculations have a hard time in extrapolating them down to such light systems as the 7r 
or the p. 

The present work opens a broad avenue of further applications, among them also the baryons. But 
thus far, it is only a first step. 
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